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extract the separated variables properly. Then by deriving the integration measure and 
the spectrum of the separated variables, we express the inner product of an on-shell and 



> 

X 

H an off-shell Bethe states in terms of a multiple contour integral involving a product of 

Baxter wave functions. Its form is reminiscent of the integral over the eigenvalues of 
a matrix model and is expected to be useful in studying the semi-classical limit of the 
product. 

^kazamaShepl . c .u-tokyo .ac.jp 
■f skomat suOhep I.e. u-t okyo .ac.jp 
^tnishimuraOhepl . c .u-tokyo .ac.jp 



Contents 

1 Introduction 2 

2 Algebraic Bethe ansatz and determinant formulas 6 

2.1 Algebraic Bethe ansatz 7 

2.2 Determinant formulas 9 

3 Separation of variables for integrable models 11 

3.1 Basic notions of the separation of variables 13 

3.2 Sklyanin's magic recipe 13 

4 Integral representation of the scalar products for XXX spin chain 16 

4.1 Construction of the separated variables 16 

4.2 Multiple integral representation for scalar products 20 

4.3 Symmetrization and simplification of the multiple integral 23 

5 Discussions 24 
A Proof of theorem 1 26 
B Relation to Izergin's determinant formula 30 

1 Introduction 

The quantum XXX Heisenberg spin chain is undoubtedly one of the most celebrated 
quantum integrable systems of fundamental interest. Besides being the prototypical model 
of the magnetic substance in the condensed matter physics, it has been under intense study 
as well from the standpoint of mathematical physics of exactly solvable models [1]. As 
such, innumerable studies have been made in the past and the basic properties of the 
model are thought to be well understood. 

However, about ten years ago this system started to receive a renewed interest from 
quite a different perspective as it made its appearance in an unexpected way: It was 
recognized as the precise mathematical structure governing the scaling properties of the 



gauge- invariant single-trace composite operators in certain sectors of the A/" = 4 SU(N) 
super Yang-Mills theory at the one- loop level [2]. In the simplest situation, where the 
composite operators are made out of two kinds of adjoint- valued complex scalars Z and 
X, Z iX) can be identified with the spin up (down) state of an individual spin forming 
an SU(2) spin chain and the dilatation operator acting on such composite operators takes 
exactly the form of the well-known XXX spin chain Hamiltonian. This identification 
allows one to compute the eigenstates and the anomalous dimensions of the composite 
operators using the techniques of the integrable models, such as the Bethe ansatz [3]. 
This technique is particularly useful when the number of magnon excitations as well as 
the number of spins become very large. In such a semi-classical limit, one can recognize 
the integrable structure quite similar to that of a classical string in a curved spacetime 
containing an AdS subspace and, as far as the spectrum of the excitation is concerned, 
this provides a structural parallel strongly suggesting the AdS/CFT correspondence [3]. 

More recently, further integrable properties of the XXX spin chain have been exploited 
beyond the spectral level to compute the correlation functions of the composite opera- 
tors of the super Yang-Mills theory. Namely, in a series of papers [4-8], the method of 
computing the three point functions ((D\{x\)02{x'2)Oj,{x-i)) has been developed, built on 
earlier works [9-11], in such a way that the three point function can be expressed in terms 
of the scalar products of the Bethe states of the spin chain Hamiltonian. For relevant 
configurations, such scalar products can be expressed in terms of the so-called Slavnov de- 
terminants [12] [13], which can be further simplified for special cases where BPS operators 
are involved. 

Now in order to study the structural similarity with the three point functions in 
the strong coupling regime represented by semi-classical string in AdS spacetime^, it 
is of importance to consider the semi-classical limit where the number of the magnon 
excitations as well as spins become very large. For the case where one of the operators 
is non-BPS and the others are BPS, such a limit has been obtained in a remarkably 
compact form [6]. More recently, the semi-classical limit of the fully non-BPS three point 
function has been worked out by Kostov [16, 17]. This was achieved again starting from 
the expression in terms of the Slavnov's determinant formula in the framework of the 
algebraic Bethe ansatz, which was quite ingenious. 

One of the strong motivations for the study performed in this paper is to construct a 
representation of the scalar products with which the semi-classical limit may be under- 



^The one corresponding to the SU(2) sector of the present interest is the string in AdSi x S^ . The 
contribution from the AdS2 part has been obtained in [14]. The contribution fr'om the S^ part, which is 
more involved, will be presented in a forthcoming paper [15]. 



stood more physically. We believe that such an understanding should be important in 
order to seek the holographic string picture in the super Yang-Mills theory. 

As has been already indicated, studies of the XXX spin 1/2 chain in the past have 

been performed predominantly in the framework of the algebraic Bethe ansatz. The 

main ingredients of this method are the elements of the monodromy matrix Q{u) = 
A{u) B{u) 



„, > T^/ \ I , which are themselves operators acting on the spin chain Hilbert space, 

with up and down spins at, say, L sites. The operators A{u) ~ D{u) satisfy the Yang- 
Baxter exchange algebra^ and from this one sees that B{u) and C{u) can be regarded 
as the "creation" and the "annihilation" operators respectively, relative to the "vacuum" 
I f^) with all the spins up. Thus one can construct the basic state with M magnon 
excitations in the form ni=i B{ui)\ '\'") and its conjugate (f-^ | nj=i C{ui). Such a state 
will be referred to as a Bethe state. It is well-known that a Bethe state becomes an 
eigenstate of the mutually commuting conserved charges, among which is the Heisenberg 
Hamiltonian of the spin chain, when the rapidities Ui of the magnons satisfy the Bethe 
equations, i.e. when they are "on-shell". Of fundamental importance in this frame work 
is the scalar product (t^ | ni=i C'(fj) nj=i B{ui)\ '\^), which can be computed using the 
Yang-Baxter algebra. When one of the set of rapidities, for instance {fj}, are on-shell, 
the scalar product can be simplified enormously and expressed as a determinant. This 
is the celebrated Slavnov determinant and practically all the calculations involving the 
scalar product have been done starting from this expression. 

In this article, however, we shall take a different route for the calculation of this 
scalar product and naturally obtain a different new representation. This alternative 
is the method of separation of variables (SoV), which was advanced substantially by 
Sklyanin [18]. The concept of SoV represents the most primitive and fundamental form of 
integrability, where one reduces the interacting many-body system to mutually decoupled 
set of dynamical systems, each with a single degree of freedom. Of course the highly 
non-trivial question is how to actually construct such separated variables {xj} and the 
corresponding canonically conjugate momenta {pi}- For the integrable systems which ad- 
mit the formulation with Lax operators, Sklyanin proposed a powerful concrete recipe for 
the construction. Relegating more detailed description to section 3, the prescription ap- 
plied to the case of XXX spin 1/2 chain says that the solutions Xi of the operator equation 
B{xi) = provide the separated coordinates, while their conjugate momenta Pi are given 
essentially by D{xi). One can indeed check that they satisfy (with appropriate ordering 
in the quantum case) the canonical Poisson-Lie commutation relations. Therefore if one 

"We shall only recall a necessary portion of this algebra later when we need them. 



can diagonalize and factorize B{u) as B{u) oc nt=i('" ~ ^i)^ with precisely L zeros, Xi 
provide a complete basis of separated coordinates. Once this is achieved, one can figure 
out the measure factor /i(xi, . . . ,xl) and compute the scalar products between various 
states in the x-representation. 

Indeed such a method has been applied successfully to some cases where the conven- 
tional algebraic Bethe ansatz is not readily applicable. One example is the non-compact 
SL(2) spin chain in the unitary representation, studied in [19]. One gratifying feature of 
this case is that in such a unitary representation the hermitian conjugate of B{u) operator 
is basically itself and hence can be easily diagonalized. The integration measure is found 
and the scalar product is thus defined in the SoV framework. Another system for which 
the SoV analysis has been performed is the SU(2) spin chain with anti-periodic boundary 

condition [20,21]. In this case, due to the insertion of the twisting matrix K = i J 

which fiips the spin at the boundary, the operator which should be diagonalized to yield 
separated variables changes from B{u) to D{u). Since this operator is hermitian and 
naturally diagonalizable the subsequent analysis a la Sklyanin is straightforward. 

Now for the more fundamental case of the SU(2) spin chain with the periodic boundary 
condition, there are two apparent obstacles in computing the scalar products using the 
Sklyanin's procedure. The first problem is that because the hermitian conjugate of B{u) 
is C{u), the basis in which B{u) is diagonal is different from the one in which C{u) is di- 
agonal. Hence the scalar product of our interest {v\u) = (f-^ | nj=i C{vi) Y[i=i B{ui)\ t^) 
cannot be easily computed in i?(M)-diagonal basis. The second problem is that B{u) 
operator as it stands is actually not a good operator in the SoV framework, since the co- 
efficient of the highest power u^ in the expansion of B{u) is proportional to S" = S^ — iS''^ 
belonging to the global SU(2), which is obviously not diagonalizable. It is perhaps for 
these reasons that this important basic model has not been treated in the SoV basis so far. 

We will solve these problems in the following manner. As for the first problem, since we 
are interested in the case where {vi} are on-shell, we may use the trick due to Kostov and 
Matsuo [22] to rewrite the scalar product into the form ~ (|^ \{S-)^-^^^^ n»=i B{vi)B{ui)\ t^), 
where only the B{u) operators appear. Then, the second problem can be solved by in- 
troducing a boundary condition changing twisting matrix K^ = i J so that the 

modified (regularized) operator B^{u) is diagonalizable. As we shall describe in detail in 
section 4, we can compute the integration measure as well as the wave functions corre- 
sponding to the general Bethe states nj=i B{wi) \ '\^) in the SoV basis. When put together 
to form the scalar product e-dependence in various quantities cancel precisely. This is as 
it should be since the original scalar product {v\u) is completely well-defined and finite. 



The original representation we obtain this way consists of contour integrals over Xj which 
surround certain poles of the integrand that depend on the index i. This can be recast 
into a more convenient form where the integration contours for all the Xj's will encircle all 
the simple poles of the integrand'^. This expression resembles the integral over the eigen- 
values of a matrix model and is expected to be useful in studying the semi-classical limit 
of the scalar product. In the Appendix A, we shall give a direct independent proof that 
our integral representation for the 2M = L case is equivalent to the Izergin's determinant 
formula for the so-called domain wall partition function of the 6- vertex model^. Using 
the Kostov-Matsuo trick and taking appropriate limits for the rapidities on both sides, 
we obtain an alternative proof that our integral representation reproduces the Slavnov's 
determinant formula for the scalar product of interest. 

The rest of this article is organized as follows. We begin in section 2 by giving a 
brief review of the framework of the algebraic Bethe ansatz where the scalar products 
of interest are defined and describe several different forms of the determinant formulas 
for them and related quantities. In section 3 the essence of the Sklyanin's method of 
separation of variables for classical and quantum integrable models will be summarized. 
With these preparations, we will derive in section 4 a new integral representation for the 
scalar product between an on-shell and an off-shell Bethe states for the XXX spin 1/2 
chain in a separated variable basis. The remaining problems to be pursued, in particular 
that of deriving the semi-classical limit from our integral representation, will be briefly 
discussed in section 5. Two appendices will be provided to give some technical details of 
the derivation. 

2 Algebraic Bethe ansatz and determinant formulas 

As a preliminary, we shall give a brief review of the algebraic Bethe ansatz for the XXX 
spin 1/2 chain and summarize the various known forms of the determinant type formulas 
for the scalar products between the Bethe states, the quantity of our prime interest. This 
will at the same time serve to fix our notations. 

■^Multiple integral formulas for the scalar product exist in the literature [23,24]. Our formula differs 
from them in form as well as in the context in which it is derived. 

^Although we will not emphasize it in this article, regarding the XXX spin 1/2 chain from the point 
of view of the 6- vertex model provides certain useful insights and there have been many interesting works 
on this topic [13,25,26]. 



2.1 Algebraic Bethe ansatz 

The basic ingredient in the framework of the algebraic Bethe ansatz is the so-called Lax 
operator acting on the product of the spin-chain Hilbert space "H and an auxiliary vector 
space. In the case of the XXX spin 1/2 chain with L sites, l-L is the tensor product of 
L copies of a two-dimensional vector space, consisting of the up-spin state | '|" ) and the 
down-spin state | ^) at each site, and the auxiliary space has the structure of C^. The 
Lax operator Ln{u) acting on the n-th site is then given by 

L„(n) ^ul + ^Y. Sy = { ^ + f « f " ) , (2.1) 



i=x,y,z ^ 



where S^^ are the local spin operators^ and u is the complex spectral parameter. We will 
impose the periodic boundary condition so that 5'^_,_^ = 5^. Going around the spin chain, 
we define the monodromy matrix Q{u) as 



Liiu-9i)-- 




M^l + m^-i 


( E SW + ^J2(^, +00 




\i=x,y,z j=l J 



^[u) = L,{u-e,)--- Ll{u ~eL)={ r>iA nL) ) (2-2) 

i^-2) . (2.3) 

Here 5* = ^^ 5^ are the total spin operators and we have introduced the inhomogeneity 
parameters 6 = {6i, . . . ,6l} at each site, which preserve the integrability. They are 
necessary for avoiding certain degeneracies in the intermediate steps and are also useful 
for other purposes^. 

Although the actions of the operators A{u) ~ D{u) on Ti are in general quite compli- 
cated and non-local, they are known to satisfy rather simple exchange relations, which we 
call Yang-Baxter algebra [1]. In particular its structure reveals that B{u) and C{u) can 
be interpreted as a "creation" and an "annihilation" operator respectively with respect to 
the pseudovacuum | f^) = | f) ® ■ ■ ■ (g) | f), in which all the spins are up. This allows one 

^ V ' 

L 

to construct the Hilbert space Ti as the Fock space spanned by the M-magnon states of 
the form \u) = B{ui)B{u2) ■ ■ ■ B{um)\ t'^), while C{y)\ t^) = 0. Mj's are the rapidities of 



u+i/2 
u-i/2' 

are generated by the operator C(f)'s as {v\ = (f-^ \C{vi)C{v2) ■ ■ -C^vm), built upon the 



the magnons, which are related to the momenta hj p = log "_^(^ . Similarly, the bra states 



'^We define S^ as 5'± = 5"^ ± iSy. 

^Although the physical meaning of the inhomogeneity parameters in the context of the three point 
functions has not been fully clarified, they are useful in generating loop corrections from the tree-level 
contributions [7,81. 



dual pseudovacuum ('|'^ | satisfying (f^ \B{u) = and (f-^ | f^) = 1- These Fock states 
will be referred to as generic Bethe states. 

Of particular importance is the transfer matrix given by T{u) = TtQ{u) = A{u) + 
D{u), which upon expanded in powers of u generates all the mutually commuting con- 
served quantities, including the Hamiltonian of the spin chain. The (dual) pseudovacuum 
is known to be the eigenstate of T[u) in the manner 

A{u)\ t^) = Qt{u)\ t") , D{u)\ t') = Q-e{u)\ t^) , (2.4) 

(t^ \Aiu) = (t" IQtiu) , (t^ \Diu) = (t^ \Q,iu) , (2.5) 

where Qe functions are defined as^ 

L L . .s 

Q,{u)^l[{u-e,), Q^{u)^l[(u-e,±'-\ . (2.6) 

Using this fact, the action of T{u) on the generic Bethe state \u) = Y[i=i B{ui)\ t'^) can 
be computed by pushing A{u) and D{u) through i?('Uj)'s using the exchange relations 
such as {u — v)A[v)B[u) = {u — v + i)B{u)A{v) — iB{v)A{u) and a similar one between 
D{v) and B{u). One then finds that \u) becomes the eigenstate of T{u) if and only if 
the following sets of equations, called the Bethe ansatz equations, for the rapidities are 
satisfied: 

This equation can also be interpreted as a periodicity condition for the phases of the 
magnon excitations as we go around the chain. When this equation is satisfied, the Bethe 
state is said to be "on-shell" (otherwise called "off-shell"). In that case, the eigenvalue 
t-a{u) of the transfer matrix T{u) is given by 

which is sometimes called the Baxter equation (2.7) for the Q-function defined as 

M 

Qu{u) = \[{u-Ui). (2.9) 



^As in these definitions, each + ( respectively — ) superscript on a function signifies that its argument 
is shifted by +| (respectively — §)• According to this convention, Q'^^{u) means Qg{u + i), etc. When 
9i = 0, the functions Qg{u) are often referred to as a{u) (for +) and d{u) (for — ). 



The equation (2.8) is equivalent to the Bethe ansatz equation through the condition that 
tu{u) has no poles despite the presence of Qu{u) in the denominator. 

The properties of the operators A{u) ~ D{u) under the global SU(2) generators S'* 
are often quite informative. For instance, from the transformation properties 

[S^Biu)] = -Biu), (2.10) 

[S+,B{u)] = Aiu)~Diu), (2.11) 

one can easily show that if the Bethe state \u) is on-shell it is the highest weight state 
with spin ^ — M. On the other hand, if it is off-shell, although having the same spin 
^ — M, it is a direct sum of states belonging to various representations and is not a highest 
weight state. 

2.2 Determinant formulas 

The main purpose of this work is to develop a method of computing the scalar product 
of the form 

M M 

{v\u) = {^'^\l[C{v,)l[B{u,)\t'^) (2.12) 

i=i j=i 

for SU(2) spin chain using the SoV formalism, leading to a new integral representation 
of such a product. Traditionally, however, the computation of such a product has been 
pursued in the framework of the algebraic Bethe ansatz reviewed in the previous sub- 
section. Although the computation is conceptually quite straightforward as one simply 
needs to move C(fj)'s all the way through B^UjYs, using the exchange algebra, and act on 
the pseudovacuum, in practice this procedure produces a multitude of terms which grow 
exponentially in the number of magnons and becomes intractable. Fortunately, in the 
case of the product between an on-shell and an off-shell Bethe states, Slavnov discovered 
a much more concise expression in the form of a determinant, which was to be called 
Slavnov's determinant formula [12]. More recently, various other types of determinant 
formulas have been developed, which are intimately related to the Slavnov's determinant. 
Since the configuration for which the Slavnov's formula is valid is precisely the one needed 
for the computation of the three point functions in the super Yang-Mills theory, which 
motivated our study, it is of interest to sketch in advance that how our new formula will 
be related, directly or indirectly, to these different variants of determinant formulas. 

As stated above, let us consider the case where either one of the set of rapidities u or 
V are on-shell. Then the original Slavnov's formula computes the scalar product {v\u) as 



a M X M determinant of the form 

rM 



UZlQti^^)Qe 




Vi-Oi- 




Uk + i 



Very recently, Kostov and Matsuo [22] showed that this expression is equivalent to an 
alternative determinantal expression of the form 

{v\u) = {-lfz'^^{z\e), z = uUv (2.14) 

where Z^^{z\d) is now a 2M x 2M determinant given by 

ii'<J^ ^ 0' \ 1=1 ^ ^ I J i<ij<2M 

They also pointed out that this equivalence is due essentially to the following equality 
valid when u or v are on-shell: 

M M 

(t^ inC'(^On^(«^)l t^) « (^' \{S-)'^-'''B{v.)Biu.)\ t^). (2.16) 

Intuitively this can be understood in the following way. Suppose the set of rapidities v 
are on-shell. Then the Bethe state nj=i -^(^i)l t^) built on the up vacuum is the highest 
weight state of global SU(2) with spin f — M. On the other hand, the state nj=i C'(fi)| i^) 
generated by the action of C{v) on the down pseudovacuum has the same eigenvalue for 
the transfer matrix T{u). Generally, an on-shell state corresponding to the same solution 
of the Bethe ansatz equations is expected to belong to the same SU(2) multiplet. Since 
Y[i=i C{vi)\ l^) is a lowest weight state with spin — ^ + M, we can make it into the highest 
weight state with spin |^ — M by the action of (5*^)^^^*^. Therefore we should have the 
equality 

M M 

l[B{v.)\ t^) oc {S^)'^"'''l[C{v.)\ ;^). (2.17) 

Taking the conjugate of this relation, we obtain (2.16). This identification will be of 
crucial importance when we develop the SoV method for the computation of the scalar 
product in section 4. 

Now it turns out that our integral formula will be more directly related to another 
variant of the determinant formula, found by Foda and Wheeler [26] . They showed that 

10 



the Kostov-Matsuo expression Z^^{z) can be identified with the so-called partial domain 
wall partition function (pDWPF) Z^^^^^{z\0), which naturally arises in the context of 
the six vertex model: 



^pDWPF(^l^) 



n 



a<l3\ 



^P)U^ 



t<j^ 



Oi 



( 



X det 



(zi-6»i+J/2)(^i-9i-J/2) 
i 

(22M-ei+j/2)(^2M-ei-i/2) 

aL-2M-\ 



{^1 



•+i/2)U-0L-i/2) 



\ 



(.ZlM~eLMl2){z2M-eL-il2) 



nL-2M~l 



. (2.1^ 



It is this expression which will be shown, in the Appendix B, to be equivalent to our 
multiple integral formula. In proving this equality, another more general determinant 
formula will be of use. It is the Izergin's L x L determinant [27] expressing the domain 
wall partition function. A slight generalization of our integral formula will be shown to be 
equal to this Izergin's large determinant and by taking the limit where L — 2M rapidities 
are sent to infinity and get decoupled, our formula and the Foda- Wheeler formula emerges 
on the respective side. As the latter can be directly shown to be equivalent to the Slavnov 
determinant [26], this proves that our formula is an alternative representation of the 
Slavnov's formula. 



3 Separation of variables for integrable models 

As reviewed in the previous section, excited states in the XXX spin chain Hi -^(^01 t^) 
are characterized as a collection of magnon excitations on top of the ground state and they 
are distinguished by a set of complex parameters called the Bethe roots, {ui}, which are 
normally interpreted as the rapidities of the magnons. Then the periodicity condition for 
such excitations leads to the Bethe equation (2.7). In this paper, however, we advocate 
an alternative view of the excited states, namely that the states are characterized by the 
nodes (zeros) of their wave functions and the Bethe roots are interpreted instead as the 
positions of the nodes. In this perspective, the Bethe equation arises as a consistency 
condition for the nodes of the wave function. 

To illustrate the basic idea, let us first discuss a simpler example, a one-dimensional 
harmonic oscillator^. As is well-known, the Schrodinger equation for the harmonic os- 



This toy model is discussed in a similar manner also in [281 
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cillator can be explicitly solved in terms of the Hermite polynomials. However, here we 
shall take a slightly different route and try to determine the spectrum without explicitly 
solving the equation. For this purpose, let us first re-express the Schrodinger equation by 
dividing both sides by the wave function ip{x): 

-M^)^^— = E. (3.1) 



2rmjj{x) dx"^ 



Then, by studying the behavior of (3.1) at large x, we conclude that il){x) should behave 
as il){x) ~ exp {—mujx'^/2h) when x is large. For excited states, ip{x) must also contain 
a polynomial prefactor, which gives rise to nodes of the wave functions. Therefore, to 
characterize ipi^) by the position of the nodes, let us write down the following ansatz for 



N 



ij{x) = l[{x-x,)e-'^^"/"'. (3.2) 

Substituting this ansatz into (3.1), we obtain the following equation 

E. (3.3) 



y^ 7 T7 7 + hu iy^ — h - I 

^ [x — Xi)[x — Xj) \^—^x — Xi 2 1 



Then from its large x behavior, the energy E is determined in terms of the number of 
nodes as i^^ = hu){N + 1/2). In addition, since the RHS of (3.3) is a constant and free 
of poles, we must demand that the residue of the poles at x = Xj on the LHS must 
vanish. This leads to a Bethe-ansatz-like equation for the positions of the nodes of the 
wave function, 

Xi = > . 3.4 

2mu ^^-^ Xi — Xj 

Although this idea of characterizing the excited states in terms of the number and 
the positions of the nodes is quite elementary and intuitive, it is technically difficult to 
apply this idea directly to the system with many degrees of freedom. However, in the 
case of the integrable models, it is often possible to decompose the system into a set of 
mutually decoupled one dimensional problems. The systematic method to carry this out 
is the method of separation of variables developed by Sklyanin, which we will explain in 
the rest of this section. By applying this method, we will see explicitly in section 4 that 
the Bethe equation for the XXX spin chain can indeed be interpreted as a consistency 
equation for the nodes of the wave function as in (3.4). 
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3.1 Basic notions of the separation of variables 

Before delving into the details of the method developed by Sklyanin, here we briefly sum- 
marize the basic notion of the separation of variables. In classical mechanics, separation 
of variable is applicable only if there are as many number of conserved charges, hi, ... , h^, 
as the dynamical variables. In such a case, by a judicious choice of canonical variables, 
it is often possible to write down a set of equations, each of which contains only one 
canonical pair {xi,Pi}: 

Wi{xi,pi;hi,. .. ,hd) = 0, i = l...d. (3.5) 

This type of equation is analogous to the expression of the energy of the harmonic oscil- 
lator, E = p^ /2m + muP'x^ jl, and one can determine the classical motion of the system 
in much the same way as in that case. 

When we consider the quantum system, the equations (3.5) are replaced by the fol- 
lowing equations for the eigenstates of the conserved charges, 

W^,(xi,p,;/ii,...,M|^) = 0, % = \...d. (3.6) 

In terms of the wave function in the coordinate representation, ^(xi, . . . ,3:^), (3.6) can 
be re-expressed as 

Wi\xi,- — \hx,...,hd\'^{xx,...,Xd) = {), i = l...d. (3.7) 

It is easy to see that (3.7) admits a completely factorized solution, \& = Yli'^ii^i)^ ^^ch 
factor of which satisfies the following one dimensional equation, 

Wi (x,, ^^; /ii, . . . , hdj ^^{Xi) = . (3.8) 

In this way, the original system with many degrees of freedom can be reduced to a set of 
mutually decoupled one dimensional systems. 

3.2 Sklyanin's magic recipe 

The most nontrivial step in the procedure above is the construction of the separated 
variables satisfying the equations of the form (3.5) or (3.6). This is indeed a difficult 
problem for interacting many-body systems. However, for the integrable models which 
can be formulated in terms of the Lax operators, Sklyanin proposed a systematic method 
for the construction, often referred to as the Sklyanin's magic recipe. In what follows, 
we sketch the essence of this recipe^ applied to systems with a 2 x 2 monodromy matrix. 



^The discussion here is basically restricted to the simplest class of the integrable models, called rational 
models. For trigonometric or elliptic models, nontrivial modification of the method is required [18]. 
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More precise analysis for the case of the XXX spin chain will be given in the next section. 

For simplicity, let us first consider the classical case. In a classically integrable system 
with a 2 X 2 monodromy matrix 

"(») - ( cS ^ ) ' l^-^' 

there is an immediate candidate for the set of equations (3.5). It is the characteristic 
equation for the monodromy matrix 

det(2-fi(x)) =0, (3.10) 

where z is the eigenvalue of the matrix VL{x). Since the expansion of VL{x) in powers of 
X yields a set of conserved charges as its coefficients, (3.10) is indeed of the form of (3.5) 
if we can somehow identify x and z with dynamical variables. The recipe proposed by 
Sklyanin is to use the solutions Xj's to the equation B{u) = as x-variables: 

B{u) = {u - xi){u - X2) ■ ■ ■ . (3.11) 

In the case of the lattice models, such as the XXX spin chain discussed in the previous 
section, B{u) is a polynomial in m, whose order basically equals the lattice size. Therefore, 
this prescription indeed provides the correct number of variables. Furthermore, owing to 
the Poisson commutativity among B{u)^s, Xj's also commute with each other and thus they 
are mutually independent separated variables. On the other hand, the 2;-variables, which 
are the eigenvalues of Q{x), are provided by the diagonal components, A{xi) or D{xi), 
since fi(xj) becomes a lower triangular matrix owing to B{xi) = 0. Then the remaining 
task is to understand the relation of A{xi) and D{xi) to the conjugate momenta p,, which 
satisfy the standard commutation relations: 

{xi,Xj} = 0, {Pi,Pj} = 0, {xi,Pj} = Sij . (3.12) 

In most cases, by explicitly computing the Poisson brackets of A{xi) and D{xi) with Xj, 
we can show that they are related to Pi roughly as 

A{xi)^e'^\ D{x,)r^e-'P\ (3.13) 

In the case of the quantum integrable models, separated variables Xj's become a set of 
commuting operators x^'s, which are characterized as the roots of the operator equation 
B{u) = 0. Just as for the classical case, the conjugate operators, e*^' and e~^^\ are given^° 



^"Note, in the quantum case, we need to consider the ordering of the operators. In the case of the XXX 
spin chain, this is exphcitly worked out in section 3. 
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essentially by A{xi) or D{xi). To derive a set of one dimensional equations of the type 
(3.7), let us consider the wave function in the Xj-basis: 

^(xi,...,Xd) = {xi,...,Xd\^), (3.14) 

where (xi, . . . ,Xd\ is an eigenstate of the operators Xj's. Now if the state |\l/) is an eigen- 
state of T{u) = A{u) + D{u), a generating function of the commuting set of Hamiltonians, 
we can compute (xi, . . . ,Xd|T(xj)|\l/) as 

{xi,...,Xd\T{xi)\^) =t{x,)^{xi,...,Xd), (3.15) 

where t{u) is the eigenvalue of T{u) for |\E'). We can evaluate the same quantity also by 
acting T(xj) to the left on (xi, . . . , Xd\- To carry this out we use the relation of T(xj) with 
the momenta pi, i.e. 

T{xi) = A(x,) + D{xi) ~ e'P^ + e-'P^ . (3. 16) 

Then we find 

(. . . , Xj, . . . \T{xi) ~ (. . . , Xi + 1, . . . I + (. . . , Xi - 1, . . . I . (3.17) 

In this way we arrive at the following equation for the wave function \&: 

t(x,)^(. . . , x„ . . .) ~ ^(. . . , Xi + 1, . . .) + ^(. . . , Xi - 1, . . .) , i = l,...,d. (3.18) 

Assuming the factorized form of the wave functions \l'(xi, . . . ,Xrf) = rijV'l^j)) "^^ can 
decompose (3.18) into a set of mutually decoupled one dimensional equations: 

t{xi)iJi{xi) ~ ij{xi + 1) + V'(xj - 1) . (3.19) 

This equation is the analogue of the Schrodinger equation for the harmonic oscillator. 
Therefore, as in that case, we can derive a consistency condition for the nodes of the wave 
function. Assuming a form of ipt as ipii^) = Yiki^ ~ ""fc) and setting Xj = Uk in (3.19), we 
obtain the algebraic relations for the positions Uk of the nodes: 

nUk — Uj + 1 , , 

— ^—— . 3.20 

. ,. Uk-Ui -1 



JT^k 



Note that this is identical with the Bethe equation. Therefore, as mentioned at the 
beginning of this section, the Bethe roots can be interpreted as the nodes of the wave 
function in this approach. In the next section, we will see that the logic outlined here is 
explicitly realized in the case of the XXX spin chain. 
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4 Integral representation of the scalar products for 
XXX spin chain 

In the preceding section, we sketched the basic idea of the method of separated variables 
for integrable models. In this section, we will apply it to the periodic SU(2) XXX spin 
chain and derive a multiple integral representation of the scalar products in the basis 
where the separated variables are diagonal. The resultant expression can be brought to a 
form which resembles the integral over the eigenvalues of a matrix model. 

4.1 Construction of the separated variables 

Recall the definition of the monodromy matrix Q{u) for the XXX spin chain with inho- 
mogeneity parameters Of. 

«(<■) = ( c[|,'| ^1",') ) s ii(« - 0,)L,(u - e,) ■■■Ldu-e,), (4.1) 

As outlined in the previous section, the separated variables for integrable models with 
a 2 X 2 monodromy matrix are usually given by the roots of the operator equation, 
B{u) = 0. However, as already pointed out in the introduction, in the case of the 
periodic SU(2) spin chain, the operator B{u) is proportional to S~ in the large u limit as 
B{u) ~ iS~u^~^ + ■ ■ ■ 5 and is not diagonalizable. This problem can be circumvented by 

introducing a twisting matrix K^ = [ J , which changes the boundary condition 

and modifies the monodromy matrix as 

a(,.)=A-.si(..)-(-(^:{:j 51:]). (4.3) 

Although such a twisting changes the dynamical properties of the spin chain, it does not 
affect the computation of the scalar products since, as we shall show explicitly later in 
this section, they can be re-expressed in terms of quantities independent of the twisting 
parameter e. After twisting, the large u behavior of B^{u) is modified to B^{u) ~ tu^ + 
i{S' —eS''^+i ^ • 9j)u^~^+- ■ ■ and B^{u) becomes diagonalizable. Then it can be factorized 
as 

L 

B,{u) = el[{u - Xi) , (4.4) 

1=1 
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where fj's are the roots of the operator equation, B^^^u) = 0. As the twisting preserves the 
algebra among the elements A{u) ~ D{u), the operators B^^s continue to commute with 
each other, namely [B^{u),B^{v)] = 0, and this implies that Xj's also mutually commute: 
[xi , Xj] = 0. These operators are the "coordinates" of the separated variables and one can 
consider their left eigenstates and right eigenstates, (xi, . . . ,xl\ and \xi, . . . ,xl), upon 
which B^{u) acts in the following way: 



{xi,...,xl\B,{u) = I e JJ(M-Xi) j {xi,...,xl\ 



(4.5) 



B,{u)\xi,...,xl) = f e JJ(n-a;j) j \xi,...,xl). (4.6) 

As explained in the Appendix A, the eigenvalue of the operator Xi takes only two values 
given hj Xi E {9i — ^ , 9i + ^Y As a consequence, the dimension of the Hilbert space 
spanned by the eigenstates of the separated variables is 2^, which precisely matches that 
of the spin chain Hilbert space. This assures the completeness of the separated variable 
basis. 

At Xi the operator B^ vanishes and the form of the monodromy matrix becomes lower 
triangular. Therefore the two eigenvalues are given by A^{xi) and D^{xi), which are 
expected to be identified as e^*^"% where Pi is the momentum operator conjugate to Xi. To 
see this more precisely, since A^{u) and D^{u) are polynomials in u with operator- valued 
coefficients, we need to specify the ordering of Xi and the coefficients, which are also 
operators in general. The ordering appropriate for the left eigenstates, to be denoted by 
: * * * -.L, turns out to be placing all the Xj's to the left of the coefficients, namely 

■.F{xiy.L = Y,^'^^n, for F(m) = ^m"F„. (4.7) 

n n 

Then the commutation relation between A^{u) and B^{u), given by {u — v)A^{v)B^{u) = 
{u — V + i)B^{u)A^{v) — iB^{v)A^{u), leads to 

:(m - Xi)A^{xi)B^{u):L = :(« - Xi + i)B^{u)A^{xi):L - i:B^{xi)A^{u):L ■ (4.8) 

Since the second term on the RHS of (4.8), containing B^{xi), vanishes, and since B^{u) 
commutes with Xi, we can simplify (4.8) to 

{u - Xi):A^{xi):LB^{u) = {u - Xi + i)B^{u):A^{xi):L , (4.9) 

where the normal-ordering is now imposed only on A^{xi). Then by acting (4.9) to the 
left eigenstate, we obtain 

L 

{u - Xi){xi, . . . , XL\:A,{xi):LB,{u) = e{u -Xi + i) JJ(n - Xk){xi, . . . , XL\:A,{xi):L ■ 



fc=i 
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Dividing both sides by (u — Xi), we see that B{u) acting on the state (xi, . . . ,XL\'-Af:{xi):L 
vanishes a.t u = Xi — i. This means that the operator ■.A^[xi):L indeed effects the shift of 
the eigenvalue of Xi by —i, namely^^ 

{...,Xi,... \:A^{xi):L (x {. . . ,Xi - i, . . .\ . (4.10) 

A similar argument for D^{u) leads to the conclusion that ■.D^{xi):L shifts the eigenvalue 

of Xi by +i, 

{...,Xi,... \:D,{xi):L oc {. . . ,Xi + i, . . .\ . (4.11) 

The constants of proportionality in (4.10) and (4.11) can be determined by the analysis 
detailed in the Appendix A. Since these results, together with the spectrum of Xi already 
quoted, are basic to the rest of the analysis, we shall display them as a theorem: 
Theorem 1: 
(i) The spectrum of Xi is given by 

i^^\e^-'-,e,+ l\ . (4.12) 

(ii) The operators ■.A^{xi):L and ■.D^{xi):L act on the left eigenstates in the following 
manner 



{...,Xi,... \:A^{xi):L = Vl + e'^Qg{xi){...,Xi-i,... |, (4.13) 

{...,Xi,... \:D,{xi):L = VlTe^Qff{xi){. . . ,Xi + i, . . .\ . (4.14) 

For the right eigenstates, an appropriate ordering prescription is to put all Xj's to the 
right of the coefficients of A^{u) and D^{u): 

■.F{xi):R = Y,Fnx'l, for F(m) = 5^m"F„. (4.15) 

n n 

Then the action of :A^{x;i):ji and :D^{xi):R on the right eigenstates are expressible as 



■.A^{xi):R\...,Xi,...) = Vl + e^Qg{xi)\...,Xi + i,...) , (4.16) 

: A(xi):fl| ...,Xi,...) = Vl + e2g+(x,)| . . . ,Xi - i, . . .) . (4.17) 

Since A^{u) and D^{u) are L-th order polynomials in u with a unit leading coefficient, the 
action of these operators at L distinct values of m, (4.13) and (4.14), completely determines 



^^For the literal identification of :A^{xi):L with e^^', it is more natural to rename Xi as —ix^. Then, 
the new Xi gets shifted by +1 and its spectrum becomes real at 6*, = 0. But we shall not do this and 
stick to the customary definition. 
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the explicit forms of the operators as follows: 

L L 



k=l k=l \jf^k ^ ^ 



L L _ 



D,{u) = W[u-Xu) + Y.[\{ —^ ) :/^e(x,):L • (4.19) 

fc=l fc=l \j^k 



•^k ^j 



They are expressible also in terms of the right-ordered operators, :A^{xi):fi and :D^{xi):ji, 
as 

A{u) = l[{u -Xk) + J2 ■■M^kY-R ( n J^] ' (4.20) 

k=l k=l \j¥=k '^ V 

D,{u) = l[{u -Xk) + J2 ■■D.XikU ( n J^] ■ (4.21) 

fc=i fc=i Vj^fc '^ v 

From (4.18) and (4.19), we can derive a difference equation for the eigenstate I?/;) 
of the (twisted) transfer matrix, T^{u) = A^{u) + D^{u). This is done by computing 
{xi,X2, ... ,xl\T^{u)\iIj) in two different ways: First by acting T^{u) on {xi,X2, . . . ,xl\ 
using (4.18) and (4.19), and second by acting it on \ip). By setting m = Xj in the resulting 
equation, we obtain the following simple equation for the wave function of the eigenstate, 

^(Xi,...,Xl) = {Xi,...,Xl\'4'). 

^'^^'^ ^{...,Xi,...)= g^(x,)^(. ..,x,-t,..:) + Qg{xi)^{. ..,Xi + t,...). (4.22) 



Here t^{u) is the eigenvalue of T^{u), i.e. T^{u)\iIj) = t^{u)\ip). Assuming a factorized 
form of the wave function, '^{xi,X2, . . . ,xl) = 4'i{xi)4'2{x2) ■ ■ .ipLi^L), (4.22) can be 
decomposed into a set of L one-dimensional equations, which can be regarded as the 
"Schrodinger equations" for the separated variables: 

^'^""'^ Mx^) = QU^i)4^r{^i) + Qgixi)^t^{x^) . (4.23) 



\/r+ 

In the e — 7- limit, the equation (4.23) for tpi apparently takes the same form as the 
Baxter equation (2.8) for the Q-function, Qu{u). However one should keep in mind that 
Qu and ipii^Xi) are conceptually quite different: While Qu is introduced as a polynomial 
with zeros at the rapidities of the magnon excitations and can be defined on the whole 
complex plane, ipi is the wave function in the separated variable basis and is defined only 
on the discrete eigenvalues Xi = 9i ±i/2. Therefore it is a priori not clear whether we can 
identify ipi with Q^. Nevertheless, as we shall later see explicitly, the factor representing 
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the wave function in the multiple integral formula is given indeed by the Q-function. 
Therefore, as far as the multiple integral formula is concerned, ipi can be identified with 
Qu and the Bethe equation can be interpreted as the consistency condition for the zeros of 
the wave function. This evidently parallels the case of the harmonic oscillator d iscussed 
in section 3. 

4.2 Multiple integral representation for scalar products 

Having constructed the separated variables, our next goal is to express the scalar product 
between an off-shell and an on-shell Bethe states given by {v\u) = (f^ | ni=i C{vi) Y[i=i B{ui) \ t'^) 
as the overlap between two wave functions of separated variables. Our basic strategy for 
deriving such a expression is to insert into the scalar product a resolution of unity in the 
SoV basis, namely 

1= Y. M=^)l=^)(^l, (4.24) 

Xi=e±i/2 

where x stands for {xi, . . . iXl} and ^{x) is the measure factor for the summation, to 
be specified later. Unfortunately, this procedure cannot be carried out straightforwardly 
because the scalar product of our interest contains the operator C{u) and its action on 
the 5-diagonal SoV basis is quite complicated. In addition, to employ the SoV basis, we 
need to introduce the twist in the boundary condition as in (4.3), which is not present in 
the original scalar product as above. 

The first problem can be circumvented by the trick due to Kostov and Matsuo [22], 
which converts C{yi) to B{yi) within the scalar product provided fj's satisfy the Bethe 
equation. Although not explicitly given in [22], one can work out the precise factors in 
the conversion formula and obtain the expression 



MM / ^ ^ M MM 



<L 



n^(^^)n^(^^)i t") = (L-2MvS ^' i(^")""'''n^(^^)n^("^)i t^)' (4-25) 



which contains only the operator B{u). This rewriting has another gratifying feature: It 
allows us to introduce the twist of the boundary condition without changing the value of 
the scalar product. This is done by replacing the second line in (4.25) with 

L MM 

{i'^\{S--eS^ + ^Y,0,f-'''l[BM)UB.inJ)\t'). (4.26) 

j=i 1=1 j=i 

Although (4.26) has apparent dependence on e as well as an extra dependence on 6'j's, 

20 



such unwanted terms actually vanisli^^ thanks to the conservation of the total spin S^ 
along the 2;-axis. 

Then, inserting a resolution of unity (4.24) into (4.26) and using the action of B^{u) 
on the SoV basis (4.5) and (4.6), we obtain the following expression^^ 

L M M 

{i' \{S- -eS^ + iY,d,f-'''X{B,{v.^\[B,{u,)\\^) 

j=l i=l j=l 

= 5^ eV(^)/L(a^)/K(^) $^2;J l[Qu{xk)Qv{xk) , (4.27) 

Xi=e±i/2 \j=l J k=l 

where fL,R are given by 

h{x) ^ {i^^ \x) , fn{x)^{x\t''). (4.28) 

Note that both the measure fi{x) and the functions fL,R{x) depend on the twist parameter 
e but the total expression (4.27) should be e-independent as argued above. 

Let us now determine ^{x) and /^ ^(a;). Note that the measure /i appears in the norm 
between the left and the right eigenstates of in the SoV basis as 

{x'\x) = ii-\x) 5^>^^^ . . . 5^>^^^ . (4.29) 

This suggests that ^{x) can be determined by computing the matrix element 

{x'\A,{u)\x) (4.30) 

in two different ways: First by acting A^{u) on the bra using (4.13) and (4.18) and second 
by acting D^{u) on the ket using (4.16) and (4.20). By equating these two expressions, 
we can easily arrive at the following recursion relation for ^{x): 



U-'^r-r^P- («i) 



/i(. ..,Xfe,...) Qt~^^ixk)W Xk-X 

The solution to this equation can be easily obtained as 

i<j k l^m ^ ' ' 

""^^To see this, it suffices to recall that B^{u) is composed of B{u) + eD{u) and S~ and that B{u) lowers 
the eigenvalue of 5'^ by 1/2 while 5'^ and D{u) leave it unchanged. 

^■^Note that the combination S^ — eS^ + ij^i ^j appears in B^{u) as Bg{u) ^ eu^ + i{S^ — eS^ + 
i ^ ■ 9j)u^~^ + . . . and its action on the SoV basis is thus given by {S" — eS^ + i ^ ■ Oj)\xi, . . . , xl) — 
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Similarly, by computing {l^ \A^{u)\x) and (a3|Ae(-u)| f^) using the formulas 

a^ \A,{u) = (i^ I {A{u) + eC{u)) = g, (n)a^ | , (4.33) 

A,{u)\ t^) = {A{u) + eC{u)) I t^) = g+HI t^) , (4.34) 

we can derive the recursion relations for fL,R{x). Solving those recursion relations, fL,R{x) 
can be determined as 

/^(a;) ocexp I -^ln(l + e2)^a;fcj , /^(a;) oc exp I ^ ln(l + e^) ^a;^ j . (4.35) 

The constants of proportionality in (4.32) and (4.35), which are left undetermined, are 
functions of the twist parameter e and the inhomogeneity parameters ^j's. These constants 
are related to the overall normalization of the scalar products and will be fixed by the 
analysis presented in the Appendix B, which compares it with the other known formula 
for the scalar product. Let us now convert the summation over the discrete spectrum of 
Xi's to contour integrals over the continuous variables Xj. This can be done by utilizing 
the pole structure in (4.32), namely 

1 _ 1 

{Xl - Omf + 1/4 ~ [xi-em+^){xi-era-^) ^^'^^^ 

Therefore, if we recall the definition of the function Qq given in (2.6), we can write the 
measure as the residue 



fi{x) = ReSz 



(4.37) 



and perform the conversion into contour integrals. Taking into account the constants of 
proportionality in (4.32) and (4.35), we finally arrive at the following multiple integral 
formula^^ for the scalar product between an off-shell Bethe-state and an on-shell Bethe 

state: 

(t^ in^(^^)n^K)i t") = '^' \l-2M)\ 



i=\ J=l 



L ^ , / L \ -^-2*^ ^ 




X 

J-J- [^ 27ri \ ^^ ' I 

k<l m 



1 ^6 v^mj'^Q Kp^m) 



In this formula, Cj denotes the integration contour which encloses 6'j±z/2 counterclockwise. 
Actually the prefactors in front of the integral are unimportant when computing physical 
observables, since they drop out upon normalizing the Bethe states. 



^''in the Appendix B, we will give a direct analytical proof of the equivalence between (4.38) and the 
known determinant formulas. 
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4.3 Symmetrization and simplification of the multiple integral 

The multiple integral formula (4.38) derived in the last subsection has one unsatisfactory 
feature: This expression becomes singular as we take the homogeneous limit, 6'j — t- 0. 
There are two sources for the singular behavior. One is that the integration contours Cj 
get pinched and collide when all the ^j's move to the origin. Another source is that at the 
same time the prefactor Y[i<ii^i~^j) ^^^^ vanish. To get around this difficulty, we wish to 
deform each integration contour into the one, to be denoted by Cau, which encloses all the 
singularities in the integrand. However, if we naively make such deformations, obviously 
we will pick up unwanted contributions coming from different integration variables encir- 
cling the poles from the same group 6i±i/2. We can avoid such contributions by inserting 
a factor of the form nfc<i(^^'^^'' — e^'^^'), which vanishes for all the undesired cases. For the 
genuine contributions for which this factor does not vanish, we must normalize properly 
to reproduce the original value of the integral. In this way, with the factor L\ coming 
from the permutation of Xj's, we obtain the following more symmetric expression for the 
scalar product: 

M M „ 



L\{L-2M)\ 




L-2M 



where the prefactor S is given by 

i<j ^ > 

Note that for this expression the prefactor H is indeed finite in the homogeneous limit. 

Although the expression above is symmetric in all the variables and hence quite useful, 
it is of interest to point out that actually we can integrate out one of the Xj's and obtain 
a slightly simpler expression containing L — 1 integration variables. To derive it, let us 
first re-express the factor nfc<i(^^™* ~ e^'^^') as a determinant of Vandermonde type: 

JJ(^2.., _ ^2..,^ ^ ^^^ (e''^^'"'^''=)i<,,fc<L • (4-41) 

k<l 

Then, by using the basic definition of the determinant, we can rewrite it into a sum over 
permutations of the form 'Y^^i^—Xy e^'^'^'^'^^^'^^^K Now note that the remaining terms in 
the integrand is completely antisymmetric with respect to permutations. Hence, all the 
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terms in the above sum contribute equally and we arrive at the following expression: 



(L - 2M)! 11/c.n 2vr^ \ tt / LV iii g^(x^)Q^(x. 



Notice that the integral is over meromorphic factors, except for exp(27r(m — l)xm)- How- 
ever for xi this factor is absent. Hence we can easily integrate out this variable by closing 
its contour at infinity. At infinity all the factors become power functions and the only 
non-vanishing integral to be performed is § dxi/{2TTixi) = 1. After this procedure, we 
may convert the factor e27r(m-i)a;m ]-,ack to the determinant and further to the original 
expression nA,<«(^^'^^'' ~ e'^'^^'). In this way we obtain the following simple formula with 
L — 1 integration variables: 



M M 



(t 



L 



nc(...)n^("i)it^)= (^_i),(;_2M), 

Note that the factor, (X]j^«)^~^*^) which was present in the previous expressions, disap- 
peared upon integration over Xi. Therefore (4.43) is structurally similar to the eigenvalue 
integral of a matrix model. Namely, Q-functions correspond to a potential term for the 
eigenvalues and nfc<i(''''fc ~ Xi)(e^'^^* — e^'"^'') can be interpreted as a modified Vander- 
monde factor. It is intriguing to note that this modified Vandermonde factor is a hybrid 
of the ordinary Vandermonde factor for the Hermitian matrix model, nfe<i(^fc ~ xi^, and 
the generalized Vandermonde factor for the unitary matrix model and the Chern-Simons 
matrix model [29], which is essentially given by nfc<i(^^'^^'° ~ e^™')^. This resemblance to 
a matrix model strongly suggests that the semi-classical limit for the scalar product can 
be analyzed by applying the method of large N expansion familiar for matrix models. It 
turns out, however, that the separated variables we deal with here and the eigenvalues of 
a matrix model behave somewhat differently and one must be very careful in adapting 
such a method. This subject will be discussed further in the final section. 

5 Discussions 

In this paper, applying the method of separation of variables a la Sklyanin, we have 
derived a multiple contour integral representation for the scalar products between an on- 
shell and an off-shell Bethe states of the XXX spin 1/2 chain with periodic boundary 
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condition. Although a compact determinant expression for such a scalar product was 
discovered already a decade ago [12] and various variants have been discussed since then, 
our novel integral formula should be useful for a number of purposes. One such area is 
the study of the semi-classical behavior of the scalar products. In particular, this is quite 
important for the comparison of the three point functions in super Yang-Mills theory 
and the corresponding dual string theory. Already results have been obtained in this 
regard [6,16,17] starting from the determinantal form of the scalar product. However, the 
procedure through which these results are obtained is ingenious but not quite systematic. 
The purpose of our present work is largely to improve on this situation. 

There are indeed apparent advantages for our formula for such a purpose. Firstly, 
almost by definition, the SoV method we employed guarantees that the factorization of 
the dependence on the basic variables occur. Such a factorization is not realized in the 
determinant formula and indeed it took some non-trivial steps to obtain such a structure 
in [16]. A further merit of the inherent factorization is that the factorized dependence 
on the inhomogeneous parameters 9i may be quite useful in applying the so-called 0- 
morphism operation [7,8], which could be a key for understanding the loop effects. 

Another advantage of our formula is that, as was shown in section 4, it can be put into a 
form reminiscent of the integral over the eigenvalues of a matrix model. An intriguing fact 
is that the measure factor of our integral formula looks like a hybrid of that of a hermitian 
matrix model and a unitary matrix model. The semi-classical limit of interest corresponds 
to the large size limit of the matrix and various techniques developed in the past may 
be utilized. However, preliminary investigation indicates that there are certain important 
differences between our integral and the matrix integral. For instance, in contrast to the 
usual matrix model eigenvalues, the separated variables x^'s do not necessarily condense. 
It appears that careful analysis for the various regions of x^'s is required in order to 
extract the semi-classical contributions systematically. Such an analysis will be presented 
elsewhere [30]. 

Apart from the semi-classical limit, our formalism should be useful for the general de- 
velopment of solvable models. An immediate application may be to the SU(3) spin chain, 
for which the SoV method is available [31,32]. It is of interest to see what modification 
from the SU(2) case would be needed to obtain a nice integral formula for the scalar 
product for such a system. 

Before ending this article, we should perhaps emphasize, alongside with Sklyanin, that 
the method of SoV, or functional Bethe ansatz as it was originally called, is in a sense 
the most fundamental and general of all the Bethe ansatz methods and hence should 
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be applicable to systems to which the other methods are not easily applicable. In this 
sense, we hope that our analysis will prove to be useful for future development of the SoV 
method. 
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A Proof of theorem 1 

In this appendix, we shall provide a proof of the theorem 1, which gives the action of the 
normal-ordered operators ■.A^{xi):L and ■.D^{xi)\L on the SoV bra (. . . , Xj, . . . |. Essentially, 
what follows is a pedagogical adaptation of the argument given in [33] . 

The proof is by mathematical induction in the number of sites L. Begin with the 
L = 1 case. The operators A{u),B{u)^C{u)^D{u) are given by-^^ 

A(u) = u-ei + iS'^ + ieS^ , B{u) = e{u-9i-iS^)+iS- , (A.l) 

C{u) = -e{u^ei+iS^)+iS^ , D{u) = u - 01 - tS'^ - teS- . (A.2) 



By solving B{xi) = for xi and substituting it into A{u) and D{u), we get 

01 + 1 



x, = e, + tS'-ie-'S-={ ^\tf ° . ) , (A.3) 



■.A{xiy.L = 2iS^ - ie-^S- + teS^ = i i ^^_^ ^^ j , (A.4) 

■.D[x^):l = -^ie + e^^)^- = -z(e + e^^) ( J q ) • (A.5) 

Since Xi is lower triangular, its eigenvalues are read off as ^i ± | and the corresponding 
normalized eigenbras (±| are given by (+| = (1,0), (— | = (l,e)/vTT~?. Then, we can 
compute the action of -.A^xi)-.]^ and ■.D{xi):l exphcitly and get 



{+\:A{xiy.L = tVl+^{-\, {-\:A{x,):l = 0, (A.6) 

{+\:D{xiy.L = 0, {-\■.D{xly.L = -^Vl + e^{+\■ (A.7) 
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We have dropped the subscript 1 for S** for simphcity. 
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This is precisely what the theorem says for L = 1. 

Next, assume that the formulas hold for up to L = A^ and consider L = N + 1 case. 
The monodromy matrix for L = A^ + 1 is given by 

^N+i = KLi ■ ■ ■ LnLm+i 1 (A. 8) 

where K is the twisting matrix given in (**). Now in order to split this into the mon- 
odromy matrix at the A^th level and the subsequent action of the Lax operator at the 
(A^ + l)th step, it is technically convenient to introduce in the final Lax operator a twisting 
matrix of a similar form, which we denote as 

K . ( i^ ;' ) . (A.9, 

with 7] being an arbitrary parameter, just like e. Therefore we write ^at+i = fijv-^jv^^i, 
where 

n^ = {KK-')L,L2---L^, KK-'=(^^^_^^ ' T "^ ) ' ^^'^^^ 

U = kL,k'\ L'^^,=KLm+i. (A.ll) 

Since the conjugation by K does not affect the structure of the algebra, we may regard 
Cli\f as the monodromy matrix for L = N for which the theorem holds. We now write the 
matrix elements of Clj^ and i^'jv+i as 

^Niu) = f 4- f-], L',,, = ( 7+^ '--^ ) , (A.12) 

\ Cat Dm I \ c^^i cLn+i J 

and compute ^Af+i = I ^ ^ ^ ^ ) • Then i?jv+i operator is given by 

\ (-'N+l Dn+i J 

Bn+i{u) = AN{u)bN+i{u) + BNiu)dN+iiu) . (A. 13) 

Let us now write the SoV basis of bras for L = A^+ 1 as {xi, . . . , Xjy; y\. By the hypothesis 
of the induction, Bn{u) is diagonal in this basis and also ^Af+i, which acts only on the 
{L + l)th site, is diagonal. Explicitly, we have 

N 

{xi,...,XN;y\BN{u) = {e - r])Y\_{u - Xi){xi, . . . ,XN;y\ , (A. 14) 

1=1 

{xi,...,XN;y\bN+i{u) =r]{u-y){xi,...,XN;y\- (A. 15) 
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We may now compute the action of Bj^j^iiu) aX u = x^ and u = y, where y is the root of 
bN+i{y) = given hj y = 6^+1 + ^5*^+1 - ir^^^Sj^^^ Since B^ixk) and &Af+i(y) vanishes 
on this state, we get 

{xi,...,X]y]y\BN^i{xk) = {xi, . . . ,XN;y\AN{xk)bN+i{xk) , (A. 16) 

{xi,...,XN;y\BN+iiy) = {xi, . . . ,XN;y\BNiy)dN+iiy) ■ (A.17) 

The RHS can be easily computed since AN{xk) shifts Xk by —i, while dN+i{y) shifts y by 
+i, with certain known factors multiplied. In this way, we obtain the formulas 

{xi, . . .,XN;y\BN+i{xk) = r]{xk - v)Qei^k){- • • ,Xfc - i, . . . ; ?/| , (A. 18) 

N 

{xi,...,XN]y\BN+i{y) = {e-r]){y-eN+i - i/2)Yl{y - Xk){- ■ ■ ,Xk, . . . ;y + i\ ■ (A.19) 

fc=i 

Having understood the action of -Bat+i at u = Xk,y on the SoV basis, we now wish 
to deduce the spectrum of Bn^i{u) using this information. Let |$) be the state which 
diagonalizes Bj^^iiu). Then by taking the inner product with the above two equations, 
we obtain 

/3(xfe)<l>(xi, ...,Xn]y)= y]{xk - r7)(5^(xfc)$(. . . , Zfc - z, . . . ; I/) , (A. 20) 

N 

/3(|/)$(xi, ...,XN]y) = {e-T]){y- 9n+i - i/2) Y[{y - Xk)<^{. . . , Xfc, . . . ; y + z) , (A.21) 

fc=i 

where $(xi, . . . , Xn', y) = {xi, ■ ■ . ,XN',y\^) and we have denoted the eigenvalue of Bn+i{u) 
by P{u). Now to simplify the analysis of the spectrum, it is convenient to extract a factor 
p{xi, . . . ,Xn', y) from $(xi, . . . ,Xn] y) in the manner 

$(a;i, ...,XN]y)= p{xu ■■■,xn; 2/)^(a;i, ...,XN;y) (A.22) 

where p{xi, . . . ,XN',y) satisfies the first order difference equations 

p{xi, ...,XN;y) = vi^k - v)pi- . . , Xfc - i, . . . ; y) , (A.23) 

N 

p{xi, . . . , xtv; 2/) = (e - r/) JJ(|/ - Xk)p{. . . ,Xk, ■ ■ .■,y + i) ■ (A.24) 



i=l 



One can easily verify that the solution to these equations is unique up to an overall 
constant. Now with such a factor removed, the reduced wave function \E' satisfies the 



^^ The uniqueness is guaranteed by the finiteness of the spectrum of Xk and y. One can construct the 
solution p by starting from the end of the spectrum. 
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equations 

/3(xfe)^(xi, ...,XN;y) = Qg{xk)'^{. . . ,Xk-i,. . ■■,y) , (A. 25) 

I3{y)^{xu ...,XN;y) = {y- On+i - i/2)^{. . . ,Xk,. . .■,y + i) . (A.26) 

It turns out that we can drastically simplify these equations by assuming the factorized 
form^^ for \E', namely 

TV 

*(xi, ...,XN;y) = xiv) n ^'^•(^'^) • (A-27) 

k=l 

The equations for \E' then get reduced to the following equations for each factor 

P{xM^) = Qt{x)U^ -i), xei9k-'-,ek+'-\ , (A.28) 

(3{x)x{x) = {y-eN+i-^/'2)x{x + i), xeleN+i-'-,0N+i+'-\ ■ (A.29) 

Note that we have used the induction hypothesis that the spectrum of each Xk is two- 
valued as above. 

The rest of the analysis is elementary. First consider the equation (A.28) and set 
X = 9k — ^- Then due to the presence of the factor Qg{x) the RHS vanishes and hence 
we must have f3{9k - ^)^k{9k - |) = 0. If ^k{9k - |) 7^ 0, then (3{9k - |) must vanish and 
6^^ — I is in the spectrum. On the other hand suppose ^fc(6'fc — |) = 0. Then ^k{9k + |) 
cannot vanish since otherwise the whole wave function vanishes. Now set x = 6'^ + | in 
(A.28). Then the RHS vanishes and so must the LHS, i.e. (3{9k + l)ik{9k + |) = 0. This 
leads to ^(^fc + |) = and hence x = 6'^ + | is in the spectrum. Similar arguments for 
(A.29) tells us that ^Af+i ± | are in the spectrum. Thus, for L = A^ + 1, we continue to 
have the same set of spectrum as stated in the theorem. 

From this analysis we learn that the finite discrete nature of the spectrum is due to two 
reasons. One is that the operators ■.A{xk)'-L and :D{xk)'-L are essentially exponentials of 
the momentum operator and hence they induce a finite shift in Xk ■ The second ingredient 
is the presence of the prefactor Qg{x). Since it vanishes at finite discrete values of x, the 
shifting must end after a finite number of steps, in the present case just one. 

What remains is the determination of the constant of proportionality in the action of 
the operators :A{xk)'-L and :D{xk)'-L- As there are only a finite number of states, such 
a constant can be adjusted rather freely by the change of the normalization of states. 
Nonetheless, there is a certain constraint coming from the following non-linear relations: 



^^This docs not miss any solution since the solution is unique. 
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N+1 



:AN+i{xky.L:DN+i{xk):L = (1 + e') H (^^ ' ^^ + V2)(£fc - ^i - 3«/2) , (A.30) 

i=l 

Af+1 

:D^+i(xfc):i:Ajv+i(xfc):i = (1 + e^) J] (3;^ - 5, - z/2)(xfc - 0, + 3^/2) . (A.31) 

These relations can be obtained in the following way. From the commutation relations 
between ■.A]\r^i{xk):L, '■DN+i{xk):L and Xk, one can show 

:AN+i{xk) --L- DN+i{xk)--L = detgflN+iixk - i/2) , (A. 32) 

■■DN+i{xk):L-AN+iixk):L = detgn7v+i(xfc + i/2) , (A. 33) 

where det^ Qn+i{u) is the so-called quantum determinant^^ , which is a central element of 
the Yang-Baxter exchange algebra. Then by using the co-multiplication rule, detq{AB) = 
detg A detq B, one can explicitly compute the RHS and obtain the relations (A.30) and 
(A.31). The constant of proportionality chosen in the theorem is compatible with these 
relations and also to the explicit equations for L = 1 case shown in (A. 6) and (A. 7) 
obtained for unit-normalized states. This completes the proof of the theorem. 

B Relation to Izergin's determinant formula 

In this appendix, we give a direct proof that a slight generalization of our new integral 
expression is equivalent to the Izergin's determinant formula [27] for the domain wall 
partition function (DWPF) which appears in the six-vertex model. From this DWPF, the 
original scalar product of our interest can be obtained by sending an appropriate subset 
of rapidities to infinity as well as requiring half of the remainder to be on-shell. 

We begin with the domain wall partition function, which is defined as follows: 

L 

ZLiw\6)^{i'^\l[Biw,)\t'^). (B.l) 

j=i 

Note that the number of B operators is equal to the number of sites L and the rapidities 

w are not restricted to an on-shell configuration. In [27], Izergin gave a determinant 

expression for this quantity, which reads 

Zl[w\U) = — '^-^= — -7- det 



U^<,i^^ - ^m - 0^) \{w. - 9, + !)(«;, - 9, - i)A<,^.<^ ■ 

(B.2) 



^^For a detailed account of the quantum determinant, we refer the reader to [1] and [33]. 
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In what follows, we shall show that this is equal to the multiple integral formula of the 
form 



a^ I n B{w,)\ t^) = ^^ \{{e, - e,m - dj + ^m - 0. 

i=l i<j 

L 



J ^) 



First, we shall transform the Izergin's formula to a form more convenient for compar- 
ison with the integral expression. By a simple decomposition, the determinant in (B.2) 
can be rewritten as a determinant of the difference of two matrices: 



Then from the definition of the determinant, we can expand the RHS of (B.5) as 

where n+ is the number of +'s in the set {cj} and the sign (—1)"+ is produced upon 
expanding the product. Now by using the definition of determinant again to re-express 
each summand back as a determinant, we obtain 

det(M^ -M+)i<„<z. = 5^(-l)"'det(M;;0i<M<i- (B.7) 

At this point, one can apply the Cauchy's determinant identity, 

to each term det(M/j)i<jj<i = det((wi — {9j — eii/2))^^)i<ij<L. Putting altogether the 
determinant (B.5) in the Izergin's formula can be expressed as 

, , , ^ , Ui<i,j<Li^i - Wj)i0j -0i- i{^j - e*)/2) 

det ' \ _ - ,j- 



(w, - 6, + 1)(^, - e, -I) J ^<^^^.<^ nt,=i(«^^ - (^3 + ^9/2) 

(B.9) 
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Substituting it into (B.2), Izergin's formula is finally transformed into the expression 



L 



Z,{w\e) = Y,{-l^\{{w.-{e,+^e,/2)) \{ ^^' ^^ ^'^'^ '^^'^\ (B.IO) 

ei=± i,j=l l<i,j<L * ^ 

which is no longer of a determinant form. 

Now we are ready to prove its equivalence to the multiple integral (B.3). This is done 
essentially by explicitly performing the contour integrals using the residue formula. By 
picking up the contributions from the zeros of the functions Qg (xj) in the denominator, 
the integral is evaluated as 



i<3 



(B.12) 



Now note the following relation, which can be checked for every pair (cj, Cj), with e, = ±1: 

(g. _ e^ + i){ei - Oj - ^ {Oj - e, - i{ti - 6,)/2) 
{e, - e, + te,){e, - Sj - ^9) i<^^ - ^j + ^(e. - e,)/2) ' 

Using this formula, the expression (B.ll) can be simplified into 

ei=± i,j=l l<i,j<L * ^ 

This is exactly the same as (B.IO), proving the assertion. 

As already stated, the original scalar product of our interest can be obtained from this 
domain wall partition function through certain manipulations. First, by sending L — n of 
the L rapidities to infinity, thereby decoupling them, one obtains the partial domain wall 
partition function with n rapidities z (2.18) in the following way: 



I j.^ ( ZL{zU{wi,...,W L-n}\0) 

[L — n)\ {wi,...,WL-n}-^oo \ iwi~^ ■ ■ ■ iw\^_^ 



^PDWPF(^I^) = ^^A_ ^^ [ ^M^ura;i,.../^W|c.J I ^3^^^^ 



If we now set n = 2M and z = uU v, where either u or v are on-shell, we recover the 
original scalar product (t"^ | ni=i C'(fi) Y[j=i -^(""i)! t'^)- O^ the other hand, if we apply 
the same manipulations to the integral formula (B.2), we obtain the multiple integral 
formula for the scalar product (4.38). This proves the equivalence of our formula and the 
determinant formula derived by Foda and Wheeler [26]. 



32 



References 

[1] V E Korepin, N M Bogoliubov and A G Izergin, "Quantum inverse scattering 
method and correlation functions," Cambridge University Press (1993) 

[2] J. A. Minahan and K. Zarembo, "The Bethe ansatz for N=4 super Yang-Mills," 
JHEP 0303, 013 (2003) hep-th/0212208. 

[3] N. Beisert et al, "Review of AdS/CFT Integrability: An Overview," arXiv : 1012 . 3982. 

[4] J. Escobedo, N. Gromov, A. Sever and P. Vieira, "Tailoring Three-Point Functions 
and Integrability," JHEP 1109, 028 (2011) arXiv: 1012.2475. 

[5] J. Escobedo, N. Gromov, A. Sever and P. Vieira, "Tailoring Three-Point Func- 
tions and Integrabihty II. Weak/strong couphng match," JHEP 1109, 029 (2011) 
arXiv: 1104.5501. 

[6] N. Gromov, A. Sever and P. Vieira, "Tailoring Three-Point Functions and Integra- 
bility HI. Classical Tunneling," JHEP 1207, 044 (2012) arXiv: 1111 .2349. 

[7] N. Gromov and P. Vieira, "Quantum Integrability for Three-Point Functions," 
arXiv: 1202.4103. 

[8] N. Gromov and P. Vieira, "Tailoring Three-Point Functions and Integrability IV. 
Theta-morphism," arXiv : 1205 . 5288. 

[9] K. Okuyama and L. S. Tseng, "Three-point functions in N = 4 SYM theory at 

one-loop," JHEP 0408, 055 (2004), hep-th/0404190. 
[10] R. Roiban and A. Volovich, "Yang-Mills correlation functions from integrable spin 

chains," JHEP 0409, 032 (2004), hep-th/0407140. 
[11] L. F. Alday, J. R. David, E. Cava and K. S. Narain, "Structure constants of planar 

N = 4 Yang Mills at one loop," JHEP 0509, 070 (2005), hep-th/0502186. 
[12] N. A. Slavnov, "Calculation of scalar products of wave functions and form factors 

in the framework of the algebraic Bethe Ansatz," Theor. Math. Phys. 79 (1989), 

502-508 
[13] O. Foda, "N=4 SYM structure constants as determinants," JHEP 1203, 096 (2012) 

arXiv: 1111.4663. 
[14] R. A. Janik and A. Wereszczynski, "Correlation functions of three heavy operators 

- the AdS contribution," JHEP 1112, 095 (2011) arXiv: 1109.6262. 

[15] Y. Kazama and S. Komatsu, to appear. 

[16] I. Kostov, "Classical Limit of the Three-Point Function of N=4 Supersymmet- 

ric Yang-Mills Theory from Integrability," Phys. Rev. Lett. 108, 261604 (2012) 

arXiv: 1203.6180. 

33 



[17] I. Kostov, "Three-point function of semiclassical states at weak coupling," J. Phys. 
A 45, 494018 (2012) arXiv: 1205.4412. 

[18] E. K. Sklyanin, "Separation of variables - new trends," Prog. Theor. Phys. Suppl. 

118, 35 (1995), solv-int/9504001. 
[19] S. E. Derkachov, G. P. Korchemsky and A. N. Manashov, "Separation of variables 

for the quantum SL(2,R) spin chain," JHEP 0307, 047 (2003) hep-th/0210216. 
[20] G. Niccoli, "Antiperiodic spin-1/2 XXZ quantum chains by separation of variables: 

Complete spectrum and form factors," Nucl. Phys. B 870, 397 (2013) arXiv : 1205 . 4537. 

[21] G. Niccoli, "Form factors and complete spectrum of XXX antiperiodic higher spin 
chains by quantum separation of variables," arXiv: 1206 . 2418. 

[22] I. Kostov and Y. Matsuo, "Inner products of Bethe states as partial domain wall 
partition functions," JHEP 1210, 168 (2012) arXiv: 1207.2562. 

[23] J. de Gier, W. Galleas and M. Sorrell, "Multiple integral formula for the off-shell 
six vertex scalar product," arXiv: 1111.3712. 

[24] W. Galleas, "Scalar product of Bethe vectors from functional equations," arXiv: 1211 .7342. 

[25] O. Foda and M. Wheeler, "Partial domain wall partition functions," JHEP 1207, 
186 (2012) arXiv: 1205.4400. 

[26] O. Foda and M. Wheeler, "Variations on Slavnov's scalar product," JHEP 1210, 
096 (2012) arXiv: 1207.6871. 

[27] A. G. Izergin, "Partition function of the six-vertex model in a finite volume," Sov. 

Phys. Dokl. 32 (1987), 878-879 
[28] N. Gromov and P. Vieira, "The AdS(5) x S**5 superstring quantum spectrum from 

the algebraic curve," Nucl. Phys. B 789, 175 (2008) hep-th/0703191. 

[29] M. Aganagic, A. Klemm, M. Marino and C. Vafa, "Matrix model as a mirror of 
Chern-Simons theory," JHEP 0402, 010 (2004) hep-th/0211098. 

[30] Y. Kazama, S. Komatsu and T. Nishimura, in progress. 

[31] E. K. Sklyanin, "Separation of variables in the quantum integrable models related 

to the Yangian Y[sl(3)]," J. Math. Sci. 80, 1861 (1996) [Zap. Nauchn. Semin. 205, 

166 (1993)] hep-th/9212076. 
[32] F. A. Smirnov, "Separation of Variables for Quantum Integrable Models Related to 

Uq{sl]sr)." Math. Phys. Odyssey 2001: Integrable Models and Beyond: in Honor of 

Barry M. McCoy 21 (2002): 455 math-ph/0109013. 
[33] E. K. Sklyanin, "Quantum inverse scattering method. Selected topics," hep-th/9211111. 



34 



